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We use a continuous-time random walk (CTRW) to model market fluctuation 
data from times when traders experience excessive losses or excessive prof¬ 
its. We analytically derive “superstatistics” that accurately model empirical 
market activity data (supplied hy Bogachev, Ludescher, Tsallis, and Bunde) 
that exhibit transition thresholds. We measure the interevent times between 
excessive losses and excessive profits, and use the mean interevent time as a 
control variable to derive a universal description of empirical data collapse. 
Our superstatistic value is a weighted sum of two components, (i) a power- 
law corrected by the lower incomplete gamma function, which asymptotically 
tends toward robustness but initially gives an exponential, and (ii) a power- 
law damped by the upper incomplete gamma function, which tends toward 
the power-law only during short interevent times. We find that the scaling 
shape exponents that drive both components subordinate themselves and a 
“superscaling” configuration emerges. We use superstatistics to describe the 
hierarchical activity when component (i) reproduces the negative feedback and 
component (ii) reproduces the stylized fact of volatility clustering. Our results 
indicate that there is a functional (but not literal) balance between excessive 
profits and excessive losses that can be described using the same body of su¬ 
perstatistics, but different calibration values and driving parameters. 
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1 Introduction 


Financial markets fluctuate as traders estimate risk levels and strive to make a profit. The 
interevent interval between times when market returns are producing excessive profits and times 
when they are producing excessive losses can be described using a continuous-time random 
walk (CTRW) formalism (see Refs, (i^) and references therein). 

Empirical market data on excessive profits and losses (5-5) define excessive profits as those 
greater than some positive fixed threshold Q and excessive losses as those below some nega¬ 
tive threshold —Q. The mean interevent time[^ between losses versus Q has been used as an 
aggregated basic variable. 

Interevent times constitute a universal stochastic measurement of market activity on time- 
scales that range from one minute to one month (5, 6). The mean interevent time can be used as 
a control variable that produces a universal description of empirical data collapse (7), i.e., the 
distribution of interevent times for a fixed mean interevent time is a universal statistical quantity 
unaffected by time scale, type of market, asset, or index. 

This distribution can be described using (i) the CTRW valley model (see Refs. (2, 4) and 
references therein), which treats time intervals as random variables, and (ii) generalized extreme 
value statistic^ for stochastic dependent basic processes {10), which are g-exponentials not ad 
hoc statistics (5, 6). Inter-event times in a multifractal structure of financial markets {11,12) 
and in the single-step memory in order-book transaction dynamics {13) are foundational in the 
analysis of double-action market activity. 


2 Principal goal 

Our principal goal is to model the empirical dataj^ associated with single-variable statistics, 
i.e., (i) the mean inter-event time period Rq between extreme (excessive) losses, defined as 
those below a negative threshold —Q, as a function of the Q{> 0) value|^and (ii) the distri¬ 
bution 'ipQ{^Qt) of inter-event times between losses Agf, previously described using ad hoc 
g-exponentials. 

Because no empirical data associated with item (i) are available, in our study of excessive 
profits we will focus on item (ii) and use the empirical data provided in Refs. (5-5). Note that 
the g-exponentials used in Refs. (5-5) cannot produce the key empirical data in item (i), and thus 
in our approach we use superstatistics. Because small losses and profits are of little concerm 
to traders, we focus on medium to high Q-values. Our goal is to provide market superstatistics 
that have universality. 

^The term ‘interevent time’ appears in the literature under such names as ‘pausing time’, ‘waiting time’, ‘in¬ 
tertransaction time’ and ‘interoccurrence time’ in the context of different versions of the continuous-time random 
walk formalism (4, 9,11-13). 

^Whether the value of losses or profits in the basic stochastic process are statistically independent is irrelevant 
because any possible correlations between them are absent in our derivations. 

^All data fits and drawings were made using Mathematica Ver. 10. 

“^For the sake of simplicity, we will treat losses as positive quantities. 
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3 Basic achievement 


We here find an analytieal elosed form of the mean interevent time period Rq between exeessive 
(extreme) losses that is greater than some threshold Q, i.e., 

fOO 

tRq^ = P{-e < -Q) = P{e > Q) = D{e)de, (1) 

where r is a time unil|^and D{e) the density of returns given by the Weibull distribution of 
extreme (or exeessive) losses (14-16), 

D(e) = 5g)’' exp (-(1)’) , e, , > 0. (2) 

Note that we eonsider random variable e to be an inerement of some underlying stoehastie 
proeess]^ Values of this random variable ean, in general, be dependent (10). Here we eonsider 
ease rj < 1 (see Table whieh means that distribution D(e) is, for e/e S> 1, a deereasing 
truneated power-law (17). 

Referenee (18) uses the Weibull distribution to deseribe the statisties of interevent times 
between subsequent transaetions for a given asset. We use the Weibull distribution and the eon- 
ditional exponential distribution of the CTRW valley model to derive superstatisties or eomplex 
statisties assoeiated with the threshold of exeessive losses or profits. 

Substituting Q into Q, we obtain 

Rq = exp ^ , (3) 


i.e.. In Rq inereases vs. the relative variable Q/e aeeording to a power-law. 

The solid eurves in Fig.[^indioate the predietions generated by ([sjj^and fit the empirieal data 
(the points are represented by different marks). This basie agreement enables us to eonstruet 
the eorresponding superstatisties and allows us to study the sueeessive empirieal data (given 
in Sees. and |^. Beeause the statistieal error is low we are able to determine t] and £ and 
derive the subsequent parameters that define the shape of the superstatisties. For example, the 
empirieal data in Fig. [^indieates that the value of the shape parameter exponent is < 1 (ef. 
Table [^, whieh for large losses or profits (i.e., e ^ e) makes the Weibull distribution Q a 
stretehed exponentially truneated power-law. In eontrast, the results presented in Ref. (5) are 
ineomplete beeause they allow no analogous eomparison with theoretieal predietions based on 
the g-exponential. 


^Later in this text we will set r = 1. 

®For the Weibull distribution the relative mean value = ^T(l/r]) and the relative variance = 


<e)2 


r(2/r)) 


- 1 


are 77 -dependent that is, they are (for fixed exponent 77 ) universal quantities. 


'Note that each curve has a slightly different multiplicative calibration parameter that dehnes its vertical shift. 
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Figure 1: Mean interevent time period Rq vs. threshold Q for four typieal elasses of quotations. 
Blaek eireles, red squares, green rhomboids, and blue triangles eoneern US/GBP exehange rate, 
S&P 500 index, IBM stoek, and WTI (the erude oil) empirieal data (from January 2000 to June 
2010), respeetively, taken from Fig. 2 in ref. (5) (plotted from the top curve down to the bottom 
one.) The solid, well fitted curves present predictions of our formula Q - values of their fitted 
parameters e and rj are shown in tab. [^(the unimportant multiplicative calibration parameter is 
not presented there). Subtle wavy deviations from these predictions are not considered in this 
work. (Empirical data were used by permission of the FPL.) 


Table 1: Values of exponent rj and quantity £ obtained from the fit of predictions of formula Q 
to the empirical data (all of them plotted in Fig. representing the exchange rate US Dollar 
against Great British Pound, the index S&P 500, the IBM stock, and crude oil (WTI). 


Index/Par. 

7 ] 

£ 

US/GBP 

0.8756U0.0156 

0.0037±0.0003 

S&P500 

0.6981 ±0.0292 

0.0035±0.0005 

IBM 

0.8246±0.0236 

0.0078±0.0007 

WTI 

0.7855±0.0182 

0.0131±0.0008 
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4 Superstatistics 


We next eonstruet an unnormalized, uneonditional distribution 'ipq^AQt) of the interevent time 
stoehastie variable, Agt, in the form of superstatistie^ based on the Weibull distribution used 
in See.|^ 

poo 

= / ijQ{AQt\e)D{e)de. (4) 

Here we assume the eonditional distribution -0^(Agile) is in the exponential formj^ 

V.g(A««k) = ^exp(-||j). (5) 

Beeause it is eonditional, the next (subsequent) loss is exaetly e, and the relaxation time is given 
by the stretehed exponential 


^q{^) = t'q(O) exp {{Bqef) 


( 6 ) 


as a straightforward extension of the exponential relaxation time used in the eanonieal version 
of the CTRW introdueed by (i, 19-21) in the eontext of photoeurrent relaxation in amorphous 
films. Here 'rQ(O) is a free (^-independent) relaxation time, and quantity Bq(> 0) is indepen¬ 
dent of variable e. Quantity Bq is a formal analog of an inverse temperature and we will later 
derive its sealing with the eontrol threshold Q value. We use the r] exponent in Q to reduee 
the number of free exponents in ([^ (Oekham’s razor prineiple) and to derive superstatisties 
ipQ{AQt) in an exaet elosed analytieal form. Note that the stoehastie dependenee of interevent 
time Aqt on loss £ assumed in Q is eonfirmed when smaller losses appear more frequently 
than larger ones. This is deseribed by definition (|^ in whieh mean time (Agt)^ = TQ{e) is a 
monotonieally inereasing funetion of e. This ereates an expanding hierarehy of interevent times 
where larger losses and profits appear less frequently than smaller ones. To make larger losses 
or profits appear more frequently than smaller ones, we ereate the opposite hierarehy of losses 
and profits using 


^q(^) = ^q(O) exp (-(Sq£)^) . 


(7) 


In this opposite ease (whieh is also analytieally solvable) we eneounter a elustering phenomenon 
in whieh shorter time intervals separate the larger values of losses/profits rather than the smaller 
ones. This eomplementary ease is briefly diseussed in See.|^ 

Substituting (|^ and Q into Q we finally derive a superstatisties in the searehed form 




= 1 Q^Q 

MQ) 


( 8 ) 


^To normalize the superstatisties given by (4 1 we divide jpglAgt) by D(s)ds = exp 


or multi¬ 


ply it by Rq. This produces conditional superstatisties limited to profits and losses no smaller then Q. 

^The exponential form of the conditional distribution (|^ assumes that the losses or profits of a fixed value e are 
statistically independent, which is generally not valid for different values of losses and profits. 
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where the sealing shape exponent 


1 

{Bqe)^ ln(rQ(£)/rQ(0))’ 

and the lower ineomplete gamma funetion 

f^qi/TQiQ) 

lEuleril + OQ, /Tq{Q)) = / exp(-|/)d|/. 

0 


(9) 


( 10 ) 


The signifieant step in the derivation of formula ([^ is the replaeement of the running vari¬ 
able e, present in the integration Q, with a new running variable y = ^Qt/Tg^e). This ehanges 
the stretehed exponential to exponential in the overall funetion under the integral in Q making 
the integration a straightforward (exaet and elosed) operation. Note that the first equality in Q 
gives a straightforward, formal generalization of the eorresponding exponent obtained within 
the eanonieal CTRW valley model (i, 19-21), where Bg is the thermodynamie f3, e is the mean 
valley depth, and the exponent value is t] = 1. 

Equation ® asymptotieally (for Agt/Tg{Q) 1) takes a power-law form 






o(Q) (Ao(/tqW))‘+“« 


Tt 


^Euler{i- + Oq) 


( 11 ) 


of the relative interevent time Agt/Tg(Q) while initially (for Agt/Tg{Q) -C 1) it takes an 
exponential form 


V’Q(AQf) 




exp 


1 + OiQ 


^gt/rgiQ) 


'^g{Q)^ + (^g V 2-f ctQ 

For ag 3> 1, Eq. ([^ reduees to the ag-independent exponential 

1 


( 12 ) 


M^gt) 


MQ) 


exp{-Agt/Tg{Q)) , 


(13) 


whieh is eonsistent with Eq. (12). 

Note that our approaeh is based solely on the relaxation time Q as a funetion of a single 
variable e. Only £ = 0, Q, is used, and parameter Q is the external eontrol threshold, i.e., by 
using (|^ and Q, we obtain 


In Rg = 


ln(rQ(g)/rQ(0)) 
lii(T-Q(e)/rQ(0)) ' 


(14) 


Thus using Q and Q we find 


'^giQ) _ 

rg{0) - • 


(15) 
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Figure 2: Collected plots of empirical distributions (colored marks) and theoretical superstatis¬ 
tics, 'ipQ{^Qt), (black solid curves), which are predictions of our formula ([^ (while the dashed 
curves were given by q-exponential shown by eq. (3) in (5)) vs. interevent time, Agf, for the 
monthly returns (a), for the relative daily price returns for sixteen typical examples of financial 
data in the period 1962-2010 (b), from minutes over the hours to daily returns for NASDAQ be¬ 
tween March 16, 2004 and June, 2006 (c), and for the detrended minute-by-minute eight most 
typical examples of financial data (d). (All empirical data were drawn from (5, 6) for permission 
of FPL and PRE.) 
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Table 2: Values of exponent aq and quantity tq(Q) obtained from the fit of formula ([^ to the 
empirieal data representing eompanies shown in Fig. [^in plots (a), (b), (e) and (d) for Rq =2, 
5, 10, 30, TO _ 


Rq 

Fig 

aq 

^q{.Q) 

Fig. 


Fig. 

«Q 

^c) 

MQ) 

Fig 

aq 

B^) 

^q{Q) 

2 

1000 

1.7699 

1000 

1.5436 

1000 

1.6129 

1000 

1.6129 

5 

3.50 

3.125 

2.30 

2.70 

3.30 

3.330 

3.60 

3.70 

10 

2.10 

5.0 

2.0 

5.0 

2.0 

4.550 

2.10 

5.0 

30 

- 

- 

1.050 

5.560 

1.0 

5.0 

1.10 

5.260 

70 

- 

- 

0.550 

4.760 

0.550 

6.670 

0.50 

5.260 


Equations (|^ and ( [T4| ) show the decisive role of the relaxation time TQ{e), but its ab initio 
derivation is difficult. Note that Fig. and Tableshow a data collapse for a given (fixed) value 
of a single control (aggregated) variable Rq. 

We analytically prove that Rq is the control variable that allows a universal form of ([^ 
that depends solely on Rq. This variable was similarly used previously in connection with the 
^-exponential (5). Using this universal form requires that we assume that the Bq in ([^ and (|^) 
depends on Q in a power-law form, or the relevant scaling relation of scaling variable Q, 


Bq = B^/^ 




£i+C 


In'^/’' R, 


Q 5 


(16) 


where prefactor B and exponent ^ are Q-independent positive control parameters. Note that the 
second equality is a scaling relation having In Rq as a scaling variable. Thus from (|^, Q, and 
(16) we obtain the superscaling of the scaling variable \nRq (or the scaling of scaling, i.e., the 


scaling of the scaling exponent). 


— = B\n^ Rq, (17) 


which we further verify by examining data, e.g., for the IBM company, which are typical of the 
empirical data used. 

From (|^, Q, and (151 we next obtain the scaling relation of the scaling variable In Rq, 


^q(Q) 

t'q(O) 


= exp 


(ln^+‘^ R, 


Sl/{l+C) 

Q 


^\n{rq{Q)/Tqm = B\n^^^ Rq. 


(18) 


Note that quantities Bq, and all depend on the single control variable lni?Q. We 

will describe and discuss the -dependence of rQ(0) below, and will use the corresponding 
empirical data to confirm all the i?Q-dependencies. 
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Table 3: Values of exponent ag and quantity tq{Q) obtained direetly from the fit of formula ([^ 
to the empirieal data represen ting IBM eompany shown in Fig. [3] for Rq=2, 5, 10, 30, and 70. 


Rq 

Q 

ag 

^qiQ) 

2 

0.0050 

1000 

1.4286 

5 

0.01389 

3.0 

3.330 

10 

0.02145 

1.90 

5.0 

30 

0.03442 

0.950 

4.550 

70 

0.04508 

0.470 

3.850 


Table 4: Universal parameter B and universal exponent C,, defining dependence of Bg vs. Q, 
obtained from the good fit of formula ( fTT] ) (solid curve) to the empirical data (black circles) 
shown in Fig. |^, for instance, for very representative IBM company. 


B 

c 

0.04798T0.0249 

2.6096±0.3478 


4.1 Empirical verification of our formulas 

Figure]^ shows the agreement between the predictions of ([^ and the empirical data for IBM for 
Rg = 2, 5, 10, 30, and 70. 

Table [3| shows the fit of quantities ag and Tg{Q). 

Figure]^a) shows the good fit of (17) (solid curve) to the data (black circles) also found in 
the third column in Tabled 

This fit allows us to determine the parameter B and exponent ( (see Table Q. 

The inset plot shows the good agreement between the ( [T^ prediction (solid curve) and the 
data (black circles). We prepared the data by putting the first equality in Q into the third 
column of Table where the e and rj of IBM is supplied in Table 

Figure l^b) shows a plot of Tg{Q) vs. Q, where Tg{Q) comes from the fourth column of 
Table 1^ The plot consists of a broken straight line or two crossing straight lines. Tableshows 


Table 5: Values of elementary quantities Bg and 'rQ(O) derived from key quantities ag and 
Tg{Q) found from the fit of formula ([^ to the empirical data representing IBM company shown 
in Fig. [^for Rg=2, 5, 10, 30, and 70 for 77=0.8246 and £=0.0078 given in tab. 


Rg 

Q 

Bq 


2 

0.0050 

0.0295 

1.4286 

5 

0.01389 

33.82951 

1.94745 

10 

0.02145 

58.86516 

1.48793 

30 

0.03442 

136.43324 

1.26807 X 10"^ 

70 

0.04508 

320.29882 

4.56615 X 10“^ 0 
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Figure 3: The superstatisties, V’Q(AQt), vs. interevent time, Agt, in log-log scale for the daily 
price returns of IBM (empty colored circles) in the period 1962-2010 for Rq=2, 5, 10, 30 and 70 
(in units of days). Black solid curves are predictions of our formula ([^ while the dashed curves 
were given by q-exponential shown by eq. (3) in (5). The dashed-dotted curves are considered 
in Sec. For Rq > 5 the power-law relaxation of 'ipQ{^Qt) is well seen for Agf > 30. The 
inset is the plot of 'ipQ{^Qt) vs. Agf in the semi-logarithmic scale for Rq=2 to clearly present 
the exponential form of the superstatisties. This exponential form was expected due to eq. (131 
as ctg is very large in this case (see tab. [^. (Empirical data were drawn from (5) for permission 


of FPL.) 


10 












Figure 4: Key dependenee of quantities: (a) l/ag and (b) tq{Q) vs. Rq obtained, for instanee, 
for the IBM eompany. Blaek eireles in main plots represent empirieal data (shown in tab. [^, 
while solid eurves are theoretieal predietions. For (a) the solid eurve was obtained by the fit of 
formula (17) to empirieal data, where the fit parameters B and (, were shown in tab. The 
indireet empirical data for the inset plot were found from eq. 0. where aq was taken from 
tab. while rj and e from tab. for the IBM company. The solid curve in this inset plot is 
the prediction of eq. ( [T^ for mentioned above parameters B and C,. For (b) the broken line or 
both solid straight lines are linear regressions (i.e. given by tq{Q) = Rq + hs, where s = 
L for the Ihs straight line and s = i? for the rhs straight line; here ro(0) = ai + hi = 1.24 as 
Rq=o = 1). Multiplicative and additive calibration parameters and hs defining both straight 
lines are shown in tab. Hence, we have an additional interpretation of tq{Q) as equal Rq up 
to some multiplicative and additive calibration parameters. The solid curve in the inset plot (i.e. 
rQ(0) vs. Rq) was obtained from formula (181, where B and Q comes from tab. while tq{Q) 
was defined by above given straight lines. 
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Table 6: Parameters of linear regressions and bg, s = L,R, defining dependenee of both 
straight lines on Rq, presented in Fig. [4|3 for the IBM eompany. 


Parameters 

L 

R 

Clg 

0.430 

-0.019 

bs 

0.80 

5.160 


the parameters of linear regressions and bg, with s = L,R, that define the dependenee of 
both straight lines on Rq. The inset plot uses ( [T5] ) and ( [T7] ) to ealeulate (i) the data points (blaek 
eireles) with tq{Q) supplied by Tableand (ii) the solid eurve, using the analytieal form of 
rqiQ)- 

Thus by analytieally and empirieally proving the i?Q-dependenee of the superstatisties 
we explain the empirieal data eollapse shown in Fig.[^ 


5 Discussion and concluding remarks 


We find an explieit elosed form of the threshold interevent time superstatisties Q that is valid 
for exeessive losses, is the foundation of the eontinuous time random walk (CTRW) formalism, 
and that is useful in the study of a double aetion market (see {13) and refs, therein). These 
superstatisties are more eredible than the g-exponential distribution that is applied ad hoc in 
this eontext in Ref. (5, 6), and they agree with the key empirieal relation between the mean 
interevent time Rq and the threshold Q (see Fig.[^. 

We model the empirieal data eollapse (ef. Fig. using superstatisties as a funetion of a 
single aggregated variable Rq and obtain, for example, the sealing shape exponent l/ag as a 


power-law funetion of In Rq and the supersealing form of (171 that is dependent upon universal 
exponent C, and prefaetor B. 

Note that @ aeeurately deseribes the empirieal statisties of exeessive profits. Here Q defines 
the threshold for exeessive profits instead of exeessive losses (see the plots in Fig. |^. We thus 
ean use the same superstatisties to demonstrate the funetional but not literal symmetry between 
exeessive losses and profits. The symmetry is not literal beeause different eontrol parameters 
B and ( and driving parameters Og and bg {s = L, R) are used. Beeause of large statistieal 
errors in the empirieal data, we eannot empirieally verify the universality of exeessive profit 
behavior. For example, for Rq = 10 exponent 1.70 < ag < 3.10 and 0.10 < tq{Q) < 0.25, 
for Rq = 30 we have 0.90 < ag < 1.50 and 0.12 < tq{Q) < 0.35, and finally for Rq = 70 we 
have 0.60 < ag < 1.10 and 0.08 < tq{Q) < 0.36, whieh exhibit ranges that are too extended. 
Note that if we substitute Tn(e) given by (|7|)— the ease opposite to that defined by — into 


Q, and use a derivation analogous to the one^rfiat produeed we obtain a result eornplemen- 
tary to @, i.e.. 


V’g(Agf) = 




X 


^q{Q) (^QVT'g(Q))^+“« 


(19) 
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Figure 5: Statistics of interevent times between arithmetic profit returns of daily closing for 
various markets (from stock exchange and forex to resource market) and time periods. All 
empirical data (discrete marks with bars) were taken from ref. (7). Solid curves are predictions 
of our formula ([^ as it can be applied both for losses and profits. Dashed curves shown, for 
instance, in plots (a), (b), and (c) are fitted by q-exponential (remaining twelve plots are very 
similarly fitted therefore, the fits are not visualized herein). However, the possible empirical 
data collapse would be incredible in this case because errors of empirical data points are too 
large. (Empirical data were drawn from (7) for permission of PRE.) 
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where 


7kzer(l + «Q, ^qV^q(Q)) = / , ^ 

“'^qV'TqCQ) 

= r£;„;er(l + Otq) — ^Euleri)- + «Q, / 'TqiQ)) i 


?/"« exp(-|/)rf|/ 


( 20 ) 


is the upper gamma funetion, whieh for Aqt / Tq{Q) 3 > 1 truneates the power-law in (19). In 
the opposite ease of AQt/TQ{Q) <C 1 we obtain 


whieh is only formally identieal to ( [TT] ). Figure shows the best predietions of the weighted 
sum of superstatisties (dashed dotted eurves) given by ([^ and ( [T9| ). Note that this eontinues to 
agree with the eorresponding empirieal data for IBM. 

We use Q, ([^, and (|^ to obtain the moment ((Agf)"*), m = 0,1, 2,... in an explieit 
elosed form, 


/o” i’Qi.^Q't)d{AQt) 
f^D(e)de 

mlfQlTQ(s)]^B(s)d£ 

!^D{e)de 

rri\ 

= RQ{TQ{Q)r- -—, m = 0,l,2,..., (22) 

1 - m/aQ 

where the first equality gives the definitior[^ Note that the moments of arbitrary order, as well 
as aq and Tq{Q), depend solely on Rq. Note also that ((Agt)™) is finite only for aq > m, 
and that otherwise it diverges. This is in eontrast to the behavior of Rq, whieh, beeause of its 
quantile (not momentum) origin, is always finite, e.g., for IBM {Aqt) is finite only for Rq < 10 
(see Table [^. We thus have two radieally different eases, finite interevent time and infinite 
interevent time {Aqt), about whieh there is mueh in the literature (see, e.g., (3, 23-25) and 
refs, therein). Infinite interevent time is partieularly interesting when it takes into eonsideration 
ergodieity breaking (26,27). 

Note that using our mieroseopie model to simulate agent behavior (28,29) gives results very 
elose to those predieted by An approaeh using agent-based modeling in this eontext was 
reeently explored by other authors (22). 

It is our hope that this work will eonstitute a strong eontribution to the researeh effort seareh- 
ing for universal properties in market behavior. 
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